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^ Abstract 

a 

d In this paper we will describe all vector-valued Siegel modular forms of degree 2 and weight 

Sym^(St)«>det*(St) with k odd. These vector-valued forms constitute a module over the ring 
of classical Siegel modular forms of degree 2 and even weight and this module turns out to 
be free. In order to find generators, we generalize certain Rankin-Cohen differential oper- 
r— I ators on triples of classical Siegel modular forms that were first considered by Ibukiyama 

and we find a Rankin-Cohen bracket on vector-valued Siegel modular forms. 



< 

1 Introduction 

^ S ^ In comparison to elliptic modular forms, Siegel modular forms and especially vector-valued 

Siegel modular forms are much less understood. Although the dimensions of the vector spaces 
of genus 2 modular forms are known due to Tsushima [ ], explicit generators are unknov^^n 
in almost all cases. For some instances however, such generators are known. These 'first few' 
I examples are due to Satoh [17] and Ibukiyama [12, 13]. 

Tsushima's dimension formula can give clues — apart from the dimensions — about the structure 
' of the modules of forms. For instance, the dimension formula can predict relations and in 

those cases the modules vdll not be free. When no relations are predicted, the modules could 
be freely generated over the ring of classical modular forms of even weight. Ibukiyama has 
. . conjectured that one of the modules he studied is in fact free, but he did not prove this [13]. In 

. ^ this paper we will prove his hypothesis and in order to do this, we develop some methods for 

^ constructing vector- valued Siegel modular forms of genus 2. Our method allows us to compute 

^ a few eigenvalues for the Hecke operators. Our results agree with calculations done by van der 

Geer based on his joint work with Faber [6, 7]. 



1.1 Siegel modular forms 

We will first introduce some notions and notation. Let V be a finite dimensional C-vector 
space, g a positive integer and let p : GL(^,C) GL(V) be a representation. A Siegel modular 
form f of weight p and genus (degree) ^ > 2 is then a V-valued holomorphic function on the 

*Korteweg- de Vries Institute for Mathematics, Universiteit van Amsterdam, Amsterdam. Current address: 
Theoretical Biology and Bioinformatics, Universiteit Utrecht, Utrecht. 
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Siegel upper half-space J^g that satisfies for every element j of the symplectic group Y g := 
Sp(2g,Z) the functional equation 

Ar-T) = pOXr,T))AT), re^g. 

Here we write 7(7, t) := ct + d for the factor of automorphy and 7 • t := (ar + 6)(ct + d)~^, where 
a,b,c,d e Ql{g,Z) and 7=|"^jerg. A Siegel modular form f has a Fourier series 

n 

where the sum is taken over the set := {n = intj) | nu e 'Z,2nij e I.,nij = nji] of all half- 
integer, symmetric g >^ g matrices n and cr(x) denotes the trace of a square matrix x. We will 
write x' for the transpose of a matrix x and when 3/ is a square matrix of appropriate size, then 
y[x] := x'yx. 

The Fourier transform a/- of a Siegel modular form f = Y,afin)q'^ defines a function a/- : Sg^ — ' y 
and if ^ > 1, then a fin) ^ =^ n>0 (the Koecher principle). Denote by Sg the subset of Sg 
of semi-positive matrices. A form f for which a fin) vanishes for non-positive n £ is called a 
cusp form. We will denote the space of modular forms of weight p by Mp{Tg) and the space of 

cusp forms by SpiTg). The general linear group GL(^,Z) ^ Tg, embedded by a |q J'-i'^' acts 
on Sg by u : n unu' and the Fourier transform af of f behaves well under this action: 

af{unu') = p{u)af{n). (1) 

From now on, we assume that g = 2 unless otherwise specified. For convenience, we often 
write square matrices |" as {a,b;c,d) and when no confusion can be possible, we will write 
(n,r/2;r/2,m) £ $2 as (n,m,r). The irreducible representations p of GL(2,C) can be characterized 
by their highest weight vector (Ai > A2) and if A2 < 0, then dimcMp(r2) = 0. Therefore we only 
have to consider 'polynomial' representations. If we write j = Ai - A2 and k = A2, then p will 
be isomorphic to Sym-'(St)i8idet^(St), where St denotes the standard representation of GL(2,C) 
and Sym-' and det* denote the j-fold symmetric product and ^-th power of the determinant 
respectively. The representation Sym-'(St) ® det^(St) is abbreviated to {j,k) and we write k to 
denote the weight (0,^). 

The structure of the ring of classical Siegel modular forms M* := ®kMk(^2) of genus 2 was 
determined by Igusa [14]. The subring := 0^=o(2)-^/fe(r2) E Af* is a polynomial ring: 

where ^4 and (pe are Eisenstein series of weight 4 and 6 (with Fourier series normalized 
at (0,0,0)) and jio and xi2 are cusp forms of weight 10 and 12 (with Fourier series normalized 
at (1,1,1)). Satoh has determined the structure of the M^-module := ®k=i(2)Mu,k)^^2) 

for ij, i) = (2,0) and Ibukiyama did the same for (j, i) = (2, 1), (4, 0), (4, 1) and (6, 0). In this paper 
we will determine the structure of Mh Our main result can be formulated as follows. 

(6,*) 

Main Theorem. The M^-module M^g ^ ^ is freely generated by seven elements Fk of weight (6, k) 
u;iiA ^£{11,13,15,17,19,21,23}. 

We shall give the elements mentioned in the above theorem explicitly, but in order to do this 
we need to introduce Rankin-Cohen differential operators on Siegel modular forms. 
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1.2 Rankin-Cohen operators 



Rankin-Cohen operators (RC -operators) send ^-tuples of classical Siegel modular forms to (pos- 
sibly vector-valued) Siegel modular forms. They are therefore useful when we want to find 
generators for modules of Siegel modular forms. RC-operators were studied in full general- 
ity by (among others) Ibukiyama, Eholzer and Choie [ , , ]. The general construction of 
RC-operators can be quite cumbersome and therefore we only explain the genus 2 case here. 

Write Hj = {p £ C[x,y] | VA e C : p(Xx,Xy) = X-^p{x,y)} for the space of homogeneous polynomials 
of degree j in two variables x and y. The representation p := Sym-'(St) is> det^(St) : GL(2,C) — 
GL{Hj) is given explicitly by: 

(p(G)-p){x,y) = det{Gf-p((x,y)G), GeGL(2,C), peHj. 

Let Rt = C[r? .|0 < i < j < 2, < s < be the polynomial ring in the 3t variables • and consider 
an element P oi Hj <s>cRt. It is convenient to write P(r];j^, ... ,r22,x,y) =: P(r-'^,. .. ,r*;t;), where 
= ('*11''*12'**12''*22) ^ symmetric 2x2 matrix and v = (x,y)'. Using this notation, we can 
give the following definition. 

Definition 1.1. An element P of Hj <^cRt is called p-homogeneous if 
PiGr^G', Gr*G'; v) = det(G)^ • P(r\ . . . , r*;G'v) 

for allGeGL(2,C). 

Ibukiyama refers to certain special elements P of the space Hj<s>cRt as 'associated poly- 
nomials', since they are 'associated' to other polynomials P called 'pluri-harmonic' polyno- 
mials [ ]. These polynomials P can be constructed from P by first choosing an element 
k = (ki,...,kt) e Z5^(j, which we will refer to as a type, and then replacing the matrices by 
with = (^? .) a 2 X 2ks matrix of indeterminates: 

P^P : Hj®cRt^Hj®cC[^'ij\0<i<2,0<j <2ks,0<s<t]. 

Note that the map P ^ P depends on the choice of the type k and that for P e Hj <s>cRt, the 
polynomial P is not necessarily pluri-harmonic. We can now give the following definition. 

Definition 1.2. An element P £ Hj ®cRt is called k-harmonic if P is harmonic in the sense 
that 

- „ a^p 

Polynomials that are p-homogeneous and ^-harmonic can be used to define RC-operators as 
shown in the following theorem and therefore we will refer to these polynomials as RC- 
polynomials . The space of RC-polynomials is denoted by ^p{k)c:Hj<s>cRt- Write t = (ti,2;z,T2) 
for an element t e ^2 and write \k \ = kg. 

Theorem 1.3 (Ibukiyama). Suppose that P e J6'p(k) with p = (],() and let fi,...,ft be 
classical Siegel modular forms on T2 of weight ki,...,kt respectively. Write dldi^ := 
(5/5Tp ^d/dz^; ^d/dz'^ yd/dT^). The H j-valued function 

Q>[P](fi, . . . , ft)(T) ■- T^r4j^ P(dldT\. . .,d/dT';v)fi(T^} ■ ■ ■ ft{T%,_,„_^,_^ , T e ^2 (2) 
(27n)'' " 
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is a Siegel modular form of weight p iS)det'^'(St) = + \k\) and genus 2. 

A more general version of Theorem 1.3 was proven by Ibukiyama [11]. Ibukiyama and oth- 
ers also give explicit examples of RC-polynomials and all RC-poljniomials for type of length 2 
were determined explicitly by Miyawaki [15]. If the length t of the type k equals 2, then a 
non-zero Siegel modular form that has been constructed using a RC-operator will only have 
weight with £ odd if a classical Siegel modular form of odd weight has been used in this 
construction [17, 12]. The cusp form ^35 £ M^^{Y2) is such a classical Siegel modular form of 
odd weight, but if we were to use it in a construction with an RC-operator, the weight of the re- 
sulting Siegel modular form will have ( > 39. The dimension of e.g. M(2,2i)(r2) equals 1, hence 
we need RC-operators with a higher type length in order to get forms of 'low' weight. This 
'trade-off causes further problems when we increase j and we will give a (partial) solution to 
this problem below by introducing a Rankin-Cohen bracket on vector-valued Siegel modular 
forms. 

Ibukiyama constructed RC-polynomials of weight (2, 1) and (4, 1) and type length i = 3 in order 
to find generators for M^^ and M^^ We will generalize these polynomials to find generators 
for M^g but our generalization can also be used to find other vector-valued Siegel modular 
forms of weight (j,/) with j > 2 and ^ > 15 odd. 

1.3 Acknowledgements 

We would like to thank Gerard van der Geer, who stimulated the author to write this paper, 
and Fabien Clery for the useful discussions and their corrections and comments. 



2 Results 



The Hilbert-Poincare series for the dimensions of M(6/)(r2) with ( odd is given by (cf [18]) 

^X:^^dimcM(6,.)(r2)-Z = (i_;,4)(i_;,6)(i_;,10)(i_^12) • (3) 

This suggests that we should look for generators for M^^^^ of weights (6,^) with ( = 
11, 13,..., 23. We will first give two methods for constructing forms of weight (6,^) with 
( = 15, 17, ... ,23 and £ = 11, 13 respectively and then we will give explicit generators. 



2.1 RC-polynomials of weight (j, 1) 

We start with an example by Ibukiyama and Eholzer [ ]. The RC-polynomials for elliptic 
modular forms were studied by Rankin and Cohen [ ] and they can be used to construct RC- 
polynomials of weight (j,0). Such a genus 1 RC-polynomial pj^k ^ C[ri,r2] can be written as 

J/2 

Pj,k^ri,r2)= )(^i+j/2 -1)^(^2 +7/2- %2_,rf-V^, j = 0(2). (4) 

j=0 

Here we use the Pochhammer s3anbol (x)n '■= x{x - l)(x -2)---{x - n + 1). The corresponding 
RC-polynomial Pj^k ^ •^(jfi)(h) is then given by 
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and it is easy to verify that Pj^k is indeed (j,0)-homogeneous and ^-harmonic. Given a polyno- 
mial p in the variables r i, . . . , r^, we then denote the map that replaces the variables by r'^iv] 
by ^: 

¥ : p(r 1, . . . , r^) " p(ri[i;], . . . , r'M) : C[r 1, . . . , r^] — 0/^-, ®ci?t • 

We will now give a construction that uses elliptic (i.e. genus ^ = 1) RC -polynomials and pro- 
duces RC-polynomials of weight (j, 1) and type k = {ki,k2,kz)- We first consider the cross prod- 
uct on the space of 2 x 2 symmetric matrices. Let J = (0, 1; -1, 0) and define for A = (a, b; b, c) 
and B = ia',b';b',c') 



AxB=AJB-BJA 



2ab'-2ba' ac'-ca' , 

, , (5) 



\ ac'-ca' 2bc'-2cb 
then iGAG')xiGBG') = detiG)-GiAxB)G' for all G e GL(2,C). 

Definition 2.1. The operator M^^ ■C[ri,r2,r3] — > ^jHj<s>cR3 is defined as follows: 

Proposition 2.2. Let p = p j^(ki+i,k2+i) be the elliptic RC -polynomial defined by equation (4), 
then M.kp is a RC-polynomial of weight {j + 2, 1) and type k = (^1,^2,^3)- 

The proof of Proposition 2.2 is elementary but very tedious. We will therefore omit the details. 
The operator M,;, can be shown to 'commute' with the Laplacian A at the cost of a shift in the 
type {ki,k2,k3) ^ {ki + l,k2 + I,k3 + 1). The above result then follows immediately. Note that 
we could also use other appropriate polynomials p e C[ri,r2,r3]. 

We now have a recipe to construct vector-valued Siegel modular forms of weight (j, with £ 
odd. First take fs e Mk^{T2) for s = 1,2,3, then choose an elliptic RC-polynomial p = 
P j-2,(.ki+i,k2+i) ^iid define P = 'M.(^ki,k2,k3)P ^ =^0,i)(^i>^2>^3)- Theorem 1.3 then tells us that 
the function ®[P](/'i,/"2, As) is an element oi M(j\^k\+i)i^2)- This would be of no use if the result- 
ing Siegel modular forms vanish identically. We can show that Qi{P]{fi,f2,fi) is non-vanishing 
by computing a non-zero Fourier coefficient. If fs = Y.n<^sin)q^ and 2{PMfi,f2,f^) = Y.nbin)q^, 
then 

b{n)= ^ P{ni,n2,n3;v)ai{ni)a2{n2)a3{n3). (6) 

(rei,re2,re3)e(SJ)^ 

Although this is a simple formula, the fact that we take a sum over all triples {ni,n2,n3) e (Sg )^ 
such that ni + n2 + n3 = n can result in a computationally difficult problem, since the number 
of partitions ni + n2 + n3 = n grows very fast as the trace a{n) grows. 

Example 2.3. The operator M.k defined in Definition 2.1 generalises the polynomials given by 
Ibukiyama [12, 13]. These polynomials are given by 

Mae^(2,ip) and Mfe((A;i + Drg -(^2 + D^i) e =^(4,i)(^) 

with k = {ki,k2,k3). Ibukiyama uses these polynomials to define Rankin-Cohen brackets on 
triples of classical Siegel modular forms and shows that the resulting vector-valued Siegel mod- 
ular forms generate the modules M^^ *) and M^^ The operator ®[My^l] can also be described 
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in terms of the cross product (5) and Satoh's RC-brackets [17]. Satoh uses the space of sym- 
metric 2x2 complex matrices as a representation space for Sym^(St) (where the GL(2,C) action 
is given hy G GAG' for A a symmetric matrix and G £ GL(2, C)) and then defines 

Let fseMk^{T2) for s = 1,2,3. Then 

F ■.= {fl,f2\^{fl,h'\eM(2 i)(r2) 

and unraveHng the definitions easily shows that F is divisible by fi in the -module M(2,*)- 
Ibukiyama's RC-brackets [/"i, /'2, /"aJ (in [ ]) are then given by [fi,f2, fa] '■= cF/fi for some non- 
zero constant c. 



2.2 A Rankin-Cohen bracket on vector-valued Siegel modular forms 

Our aim is to find all Siegel modular forms of weight (6, £) with £ odd. Formula (3) shows that 
dimcM(6,ii)(r2) = 1 and the lowest £ for which we can use the operator Mk to construct a non- 
zero form of weight (6,^) equals 15. Note that by taking k = (4,4,4) and p = P4,(5,5), we can get 
a form 2i[M.kp]{(p4,(p4,(p4) £ M(6,i3)(r2), but unfortunately this form vanishes identically. This 
means that the above described construction with M.^ is not sufficient for our purpose, unless 
we would be able to divide by a classical Siegel modular form, but this appears to be quite 
difficult without prior knowledge of the quotient. Ibukiyama encountered a similar problem 
when he determined all modular forms of weight (6,^) with £ even and he solved this by using 
a theta series ©s of weight (6,8) (cf [13, 9]) and a Klingen-Eisenstein series Eq of weight (6,6) 
(cf [ ]). These forms can not be constructed directly by means of RC-operators. However, as 
we will point out below, we can use RC-operators to compute their Fourier coefficients. In this 
paper we have scaled ©g such that the Fourier coefficient of ©g at (1, 1, 1) equals x'^y^ + 2x^y^ + 
x^y"^. The form Eq is scaled such that its Fourier coefficient at (1, 0, 0) equals x^. 

The forms ©s and Eq can be used to construct forms of weight (6,13) and (6,11). In order to 
see how this can be done, we must first give a new interpretation of the operators Q[M.kp], 
where p is an elliptic RC-polynomial. Choose p = p j-2,(ki+i,k2+i) and let fs £ Mk^{T2) for s = 
1,2,3. Also let q = pjxkiM) and define F = ®[^g](/"i,/"2) £ M(^j^ki+k2)(^2)- We can re-write G := 
Q>[Mkp](fi,f2,h) in such a way that 

G = c-{F,h} 

for some bilinear form {•,•} on M(j^ki+k2)^^'i) x M.k^iT2) and a constant c £ C*. We will give 
the exact description of {•,•} below, but let us first state the advantage of this effort. We can 
replace the modular form F that was defined via a RC-operator by any other modular form of 
weight (j,k\ + k2). So for instance, we can take F = Eq and fs = (p4 and if we then apply {•,•}, 
we get 

{£;6,^4}eM(6,ii)(r2). 

We will now give an explicit formula for {•,•}. 

Definition 2.4. Suppose that F £ M(j^k)i^2) and (p £ Me{T2)- Define the determinant W(r) by 



W(r) ■-- 



rii ri2 r22 

9 9 

y -xy X 

dx^ dxdy dy^ 



■.Hj^Hj^cRi- 
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We now define the brackets {•,•} as follows. 



(ph^i (^^ ^^^^ " ^^^^ - ^'^^^^)^ 

By our considerations above, we then have the following result. 

Proposition 2.5. Take F and (p as in Definition 2.4. The form {F,(p} is then a Siegel modular 
form of weight ij,k + ( + 1). 

In order to be able to work with the forms {i^, we need to know how to compute Fourier coef- 
ficients. We can easily find a formula similar to (6). Write ip = Y.n>oO-MQ'^ a^id F = T.n>o b(n)q"' 
and let {F,(p} = T.n>oc(n)q"' , then we get 

0-l)c(n)= X ik+jl2-l)a{ni)W{ni)b(n2)-ea{ni)W(n2)b(n2). (7) 

(ni,ra2)£(SJ)2 
ni+re2=n 

Remark. A modular form that is defined using the brackets {•,•} will always be a cusp form. 
Suppose that we use F e M{j^k)i^2) and ip e M(iT2) in order to construct [F, ip] e M(j^k+e+i)i^2), 
then either ^ + ^ + 1 is odd and hence [F, (p\ must be a cusp form, or one of the integers £ or k 
is odd, implying that F or ^ is a cusp form. 

We can also see directly from Formula (7) that the brackets map to S{T2) by considering the 
Fourier coefficients at singular indices n. Suppose that n = ni + n2 where n,ni, n2 e and n is 
singular, then we can find a ue SL(2,Z) such that unu' = (v,0,0), uniu' = (vi,0, 0) and un2u' = 
(v2>0,0). Hence, by Formula (1) we can assume without loss of generality that n is of the form 
n = (v,0,0). The Fourier coefficients b{n) oi F for n = (v,0,0) always have the form a{v)-x-^ for 
some constant a(v). Therefore, we get for s = 1,2 that 



W{ns)b{n2} = Vs 



-xy y? 



a(v2)^ = 0. 



dxdy dy^ 

This shows that c{n) in Formula (7) vanishes for singular n. 



2.3 Generators for M}^ . 

(6,*) 

As promised, we will now give explicit generators for M^g We first use the brackets {•,•} to 
define forms of weight (6, 11) and (6, 13): 

Fii := {Ee,ip4}/1152 , 13 := {©s, <P4}/4 

and for the remainder of the generators, we use the construction with M.^. Choose the following 
elliptic RC -polynomials pi and types ki: 



Pi 


polynomial 


type ki 


Pl5iri,r2) = 


= (Brj- 


14rir2 + 7r^)/160 


(5,4,5) 


^17(^1,^2) = 


-- i4rl- 


8rir2 + 3r|)/192 


(5,6,5) 


Pl9(ri,r2) = 


-- {22r\ 


-24nr2 + 5r|)/1920 


(4,10,4) 


P2l(ri,r2) = 


-- (22r\ 


- 24^1^2 + 5r|)/2880 


(4,10,6) 


P23(ri,r2) = 


-- (ISrf 


-14rir2+3r|)/16 


(5,12,5) 
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k 


A(2) on A^(6,«(r2) 


A(2) on S(6,«(r2) 


6 


-24 -(1 + 24) 




8 







10 


216-(l + 2^) 


1680 


11 




-11616 


12 


-528 •(1 + 2^°) 


Z2-22368Z + 57231360 


13 




-24000 


15 




+ 68256X + 593510400 


17 




+ 363264^2 + 136028160Z - 4603543289856000 


19 




+ 1202400^3 - 1311202861056^2 
-179858880190218240X - 1566691549034368204800 



Table 1: Eigenvalues of the Hecke operator r(2) on M(e^k)i^2) for some values of k. If a polynomial in X 
is given, the eigenvalues A(2) are the roots of this polynomial. The space A^(6,^;)(r2) is the orthogonal 
complement of <S(6,^;)(r2) with respect to the Petersson product. 

We then v^rrite Pi = M^^.p; for i = 15, 17, ... ,23 and define modular forms Fi as follow^s. 



Fi 


modular form 


weight 


Fi5 -- 




(6,15) 


Fn -- 


= si[PnMx5,na5) 


(6,17) 


Fi9 -- 


= ^[Pl9M(P4,XlO,(P4) 


(6,19) 


F21 = 


= ^[P2l'\i(P4,XlO,(P6) 


(6,21) 


F23 = 


= ®[-P23](l5,I12,l5) 


(6,23) 



The form X5 denotes the square root of xio in the ring of holomorphic functions on (see 
e.g. [8]). Our main result can now be formulated as follow^s: 

Main Theorem, The M^-module M^g is free and can be written as the following direct sum: 

iEl 

where Fi are defined above and the direct sum is taken over the set I = {11, 13, 15, 17, 19,21,23}. 
2.4 Eigenvalues of the Hecke operators 

Since we now know all modular forms of weight (6,^) vdth ^ e Z, we can calculate the eigen- 
values of the Hecke operators Tip) (cf [2, 1]). We did this for p = 2,3 and some k (Table 1 
and 2). 

At our request, G. van der Geer computed some of these eigenvalues using a completely 
independent method that is based on counting points on hyperelliptic curves over finite 
fields [9, 6, 7]. The values listed here agree with these. Also note that the characteristic 
poljmomials of T(2) and T(3) on S(^e,k)(^2) with ^ = 12 and 15 have the following discriminant: 



k 


A(det(T(2)- 


X)) 


A(det(T(3)-Z)) 


12 


2103272301 




2i43672;L32601 


15 


2^03^29 -83 


103 


2^23^29 -53283 -103 
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k 


A(3) on S(e,ki^2) 


8 


-27000 


10 


-6120 


11 


-106488 


12 


+ 335664Z - 14832719455680 


13 


-8505000 


15 


Z2+228022128X + 8319716602228800 


17 


+ 1086146712^2 - 341960280255362880Z - 188775313801934579676864000 



Table 2: Eigenvalues of the Hecke operator r(3) on M(Q^k)(T2) for some values of k. 



This shows for p = 2,3 that the eigenvalues of T{p) on S(6,i2)(r2) and S(6,i5)(r2) are elements 
of the same quadratic number field (Q'(\/601) and Q(\/29-83- 103) respectively. We also verified 
that the characteristic polynomials of T(2) and T(3) on S(6,i7)(r2) define the same number field. 

3 Proof of the main theorem 

We w^ill only give a sketch of the full proof since it involves many elementary computations. 
Let U denote the Hodge bundle corresponding to the factor of automorphy j. The forms Fi 
defined above are sections of Sym^(i7)<8iL®^', where L denotes the line bundle det(f/). In order 
to show that the forms Ft axe independent over M^, we have to show that xi40 := ■'''ii Ai*'i3 A 
... AF23 eL®^^" is non-vanishing. The form is an element of Mi4o(r2) and we can compute 
the Fourier coefficients of this form. We will therefore prove the main theorem by computing a 
non-zero Fourier coefficient of X140. Write juo = T.nCin)q" and Fi =Y.n<^i^n)q^ . The following 
formula holds for the Fourier coefficients c(n): 

c(n)= det(aii(ni),ai3(?i2),...,a23('i7))- (8) 

As mentioned above, this can be hard to compute if cr(n) is large. In order to find a non-zero 
Fourier coefficient c{n) of X140, we need at least o{n) > 14 since the forms Fi are cusp forms. 
Indeed, under the Siegel operator the forms Fi map to elliptic modular forms of odd weight. 
These forms all vanish. 

Using algorithms provided by Resnikoff and Saldana [16], we calculated Fourier coefficients 
of the classical Siegel modular forms (p4,(p6,Xio and Xi2- The Fourier coefficients of js can be 
computed using jg = xw-^ We then were able to compute Fourier coefficients of the forms Fi. 

In order to find Fourier coefficients of Fn and F13, we first had to compute Fourier coeffi- 
cients of and ©8. The modular form (P4EQ is an element of M(6,io)(r2) and this space has 
dimension 2. Using a RC-polynomial of weight (6,0) and the modular forms (p^ and (pe, we 
can find a modular forms Fio in the space M(6,io)(r2) (cf Ibukiyama [ ]) and this form is not 
an eigenform. Hence, we can find Fourier coefficients of i*'io and T(2)i*'io. For some a,/i e Q, 
we must have aFio + fiTi2)FiQ = cpiEe. The form is an eigenform and by a theorem due to 
Arakawa [2] we know the corresponding eigenvalue of T(2) (cf Table 1). This gives a relation 

-"^ While computing the Fourier coefficients of X5 we encountered an error in Table IV of [ ]. The coefficients at 
calculation classes (3,3,3) and (2,6,0) should have opposite sign. 
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i 


11 


13 


15 


17 


19 


21 


23 


n 


(1,1,0) 


(1,1,1) 


(2,1,0) 


(2,1,0) 


(2,1,1) 


(2,1,1) 


(2,2,1) 


ai(n) 














1 


-5 


3 




-20 


-2 


312 





14 


-10 


-37 







-5 








36 


-6 


-50 










180 


-300 


24 


-24 










5 











-30 


50 




20 


2 


-102 


354 





-12 


37 






















-3 



Table 3: A few Fourier coefficients of the forms Fi. The Fourier coefficients are written as column 
vectors. The determinant of the above matrix occurs in the sum (8) for c(12,8,4). This determinant 
equals 21^3^5^11. 

for the Fourier coefficients of E% and using this relation, we were able to find a and (i. This 
also allowed us to compute the Fourier coefficients oiE%. 

A similar method can be used to compute Fourier coefficients of Qg- Again using a RC-operator 
of weight (6,2) and the modular forms (p^ and (pQ, we can find a modular form F\2 in S(6,i2)(r2) 
(cf Ibukiyama [13]). This form is again not an eigenform. The form (piQ% £ 'S(6,i2)(r2) must 
be a linear combination aFi2 + pT{2)Fi2 and since Ibukiyama has computed a few Fourier 
coefficients of ©8 (cf [ ]), we were able to determine a and p. 

We then used formulas (6) and (7) to compute Fourier coefficients of the forms Fi. A few exam- 
ples are given in Table 3. 

We wrote a script to compute the Fourier coefficient of xuo at n = (12,8,4) and found that 
c(12,8,4) = -2^^3'^5^ # 0. We checked our computations by also computing the Fourier coeffi- 
cient at (12,8,-4) and found that c(12,8,-4) = c(12,8,4) which is in line with equation (1). This 
proves our main result. 
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